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Abstract —  In  this  research,  we  have  developed  an  algorithm  to 
reduce  the  residual  artifacts  of  the  background  clutter  (that  is, 
stationary  targets)  that  appear  in  the  MTI  imagery  that  are 
generated  by  Global  Signal  Subspace  Difference  (GSSD)  of  the 
monostatic  and  bistatic  images  of  an  along-track  monopulse 
synthetic  aperture  radar  (SAR)  data.  We  have  also  established 
the  theoretical  foundation  for  estimating  the  motion  track  and 
parameters  of  the  detected  moving  targets.  We  will  show  the 
results  of  these  algorithms  on  measured  SAR  data. 

I.  Introduction 

In  recent  years,  the  United  States  Air  Force  Research 
Laboratory  at  Wright  Patterson  Air  Force  base  has  conducted 
nonlinear  SAR  data  collections  with  an  X-band  radar  system 
[1][5].  This  platform  is  capable  of  making  bistatic 
measurements  in  the  along-track  domain  as  well  as 
monostatic  measurements.  The  data  were  collected  over  an 
urban  area  that  contained  high  stationary  clutter  (buildings, 
vegetation,  etc.).  The  radiated  scene  also  included  various 
roads  with  different  types  of  moving  vehicles. 

Due  to  the  heavy  clutter  in  the  imaging  scene,  the 
signatures  of  most  the  moving  targets  fall  below  the  clutter 
signature.  We  have  developed  a  nonlinear  subaperture-based 
coherent  processing  of  dual  receiver  channels  of  the  Gotcha 
platform  to  detect  the  moving  targets  [4].  This  approach 
utilizes  the  2D  adaptive  processing  that  we  proposed  in  [2]- 
[4]  to  blindly  calibrate  the  two  receiver  channels  of  an  along- 
track  monopulse  SAR  system.  Furthermore,  a  pre-processing 
of  the  nonlinear  subaperture  data  ensures  the  two  channels  are 
coherently  calibrated  with  respect  to  the  variations  of  the 
flight  path.  The  current  paper  is  concerned  with  calibration 
of  the  two  along-track  monopulse  SAR  images  via  a  global 
spatially -varying  2D  adaptive  filter  that  provides  more 
accurate  information  not  only  for  ground  moving  target 
indication  (GMTI)  but  also  tracking  data  that  could  be 
exploited  to  estimate  the  motion  parameters  of  a  moving 
target.  The  basic  principle  behind  extracting  moving  target 


indication  (MTI)  from  two  monostatic  and  bistatic  along- 
track  receiver  channels  of  a  SAR  platform  is  that  after 
compensating/calibrating  for  known  deterministic  differences 
of  the  two  channels  the  SAR  MTI  image  can  be  constructed 
via  the  following: 

fun  L>')  =  f m  (x,  y)  -  fb  (x,  y ) , 
where  fmix,y)  and  fb{x,y)  are,  respectively  the 

monostatic  and  deterministically-calibrated  mono  static  and 
bistatic  SAR  images.  In  practice,  due  to  unknown  variations 
of  the  electronics,  antennas,  etc.  of  the  two  receiver  channels, 
there  are  unknown  phase  and  gain  variations  in  both  range 
and  Doppler  domains  that  are  unknown  to  the  user.  The 
simplest  way  to  model  this  is  to  assume  that  these  variations 
are  invariant  in  range  and  Doppler.  In  that  case,  under  the 
null  hypothesis  that  is  there  is  no  moving  target,  the 
monostatic  and  bistatic  images  are  related  via  the  following: 

fm{x,y)  =  fb{x,y)®h{x,y) 

=  j/*b ~u,y  —  v)  h(u,v)dudv 

where  ®  represents  two-dimensional  convolution, 
and  A(x,  y)  is  an  unknown  two-dimensional  filter.  Under 
the  null  hypothesis,  this  filter  can  be  determined  using  the 
least  mean  square  (LMS)  algorithm;  this  approach  is  called 
adaptive  filtering.  A  practical  implementation  of  this  method 
for  the  two-dimensional  problems  was  described  in  our 
previous  work  [2] -[4],  and  was  referred  to  as  Signal  Subspace 
Processing  (SSP). 

A  more  realistic  miscalibration  model  for  the  two 
receiver  channels  is  based  on  the  fact  that  the  filter  is 
spatially-varying.  In  this  case,  the  relationship  between  the 
monostatic  and  bistatic  images  can  be  expressed  via  the 
following: 
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fm(x>y)  =  \fb(x-u>y-v)  hxy{u,v)du  dv 

where  in  this  model  the  filter  hxy  (u,v)  varies  with  the 

spatial  coordinates,  that  is,  (x,  y) .  While  the  above  model  is 

a  more  suitable  one,  however,  it  is  computationally 
prohibitive  to  implement  the  LMS  or  SSP  method  for  this 
scenario. 

A  practical  alternative  is  to  assume  that  the  filter  is 
approximately  spatially-invariant  within  a  small  area  in  the 
spatial  domain.  In  this  case,  we  can  divide  the  SAR  scene 
into  subpatches  within  which  the  filter  can  be  approximated 
to  be  spatially-invariant.  The  resultant  model  is  given  by: 

fmt  (x>  y)  =  fu  (x>  y ) 0  K  (x>  y) 

=  J/m(x  — y~v)  he(u,v)dudv 

where  £  represent  an  index  for  the  subpatches. 

In  the  approach  that  we  call  Local  Signal  Subspace 
Processing  (LSSP),  the  LMS/SSP  method  is  used  to  estimate 

the  local  unknown  calibration  filter  h£(x,y).  After  this 

filter  is  estimated  for  each  subpatch,  an  approach  that  we  call 
Global  Signal  Subspace  Processing  (GSSP)  is  used  to 
estimate  the  original  spatially- varying  filter  hxy  {u,v). 

II.  Mathematical  Prelimenaries 

The  basic  foundation  of  our  work  is  the  use  of 
subaperture  processing  to  not  only  detect  but  also  geolocate 
and  track  moving  targets.  The  simple  and  most  prominent 
rational  for  exploiting  subaperture  processing  is  that  the 
signature  of  a  moving  target  appears  fairly  localized  in 
subaperture  SAR  imagery  though  it  is  slightly  smeared  (Note 
that  the  signature  of  a  moving  target  in  a  full-aperture  image 
appears  as  a  severely  smeared  structure  in  the  formed 
image.).  We  should,  however,  point  out  that  the  subaperture 
image  of  a  moving  target  may  not  be  visible  in  the  formed 
image.  In  fact,  the  coherent  processing  of  the  monostatic  and 
bistatic  images  is  the  key  in  suppressing  the  background 
clutter  and  revealing  the  signature  of  the  moving  targets  in 
the  imaging  scene. 

We  should  also  point  out  that  subaperture  processing  has 
also  merits  in  the  case  of  SAR  data  that  are  collected  over  a 
nonlinear  flight  path.  The  reason  for  this  is  not  just  for 
imaging  purposes  for  wavefront  reconstruction.  It  turns  out 
that  a  nonlinear  subapertures  that  result  in  interrogating 
(imaging)  the  scene  at  different  slant  planes  actually  adds  to 
the  information  base  that  is  exploited  to  estimate  the 
parameters  of  a  moving  target.  This  will  be  discussed  later. 

This  section  is  intended  to  outline  the  basic  equations 
that  model  the  signature  of  a  moving  target  in  a  subaperture. 
This  model  is  the  foundation  of  our  approach  to  geolocate 
and  track  moving  targets.  We  denote  the  slow- time  domain 


with  T  .  We  divided  the  SAR  data  into  subapertures  that 
overlap.  The  size  of  a  subaperture  is  chosen  such  that  the 
basic  assumption  that  the  subaperture  image  of  a  moving 
target  is  relatively  localized.  This,  clearly,  depends  on  the 
SAR  system  parameters  (frequency  band,  platform  speed, 
etc.)  and  the  anticipated  maximum  speed  of  the  ground 
moving  targets.  In  the  case  of  the  Gotcha  platform,  a  suitable 
subaperture  size  is  about  1,204  PRIs.  In  our  implementation, 
we  use  512  PRIs  for  the  overlap  region  between  two  adjacent 
subapertures. 

Consider  the  £  -th  subaperture  of  the  flight  path.  In  the 
subaperture  wavefront  processing,  a  spatially-varying  motion 
compensation  is  used  to  convert  the  measured  data  over  the 
nonlinear  path  every  subaperture  into  a  conventional  linear 
SAR  with  a  constant  platform  velocity  for  that  subaperture; 
we  call  this  the  synthesized  linear  SAR  subaperture.  (We 
should  emphasize  that  this  approach  motion  compensates  the 
measured  SAR  data  for  every  individual  point  in  the  3D 
spatial  domain;  that  is,  it  is  as  accurate  as  the  backprojection 
imaging.)  We  denote  the  radar  coordinates  at  the  midpoint  of 

the  resultant  linear  SAR  for  the  £  -th  subaperture  by 

v(t)  y(0  y(£)  ) 

We  also  denote  the  constant  velocity  of  the  synthesized  linear 
SAR  subaperture  by 

(v(£)  V(£)  V(0  ) 

V  xradar  ’  yradar  ’  zradar  / 

Note  that  the  synthesized  platform  velocity  does  vary  from 
one  subaperture  to  another. 

Next,  we  identify  the  slant-plane  where  the  image  is 
formed  for  the  £  -th  subaperture.  We  denote  the  3D  spatial 
domain  by  (x,y,z).  (This  is  equivalent  to  what  is  called 
PCS  coordinate  system  by  General  Dynamics.)  The 
wavefront  processor  forms  a  slant-plane  image  with  motion 
compensation  that  incorporates  variations  in  terrain  (using 
Digital  Elevation  Map,  DEM).  The  final  output  of  the 

wavefront  processor  is  a  2D  image  on  the  (x,jp)  plane  that 
is  interpolated  from  the  slant-plane  SAR  image.  For 
convenience,  we  refer  to  the  2D  ground  (x,y)  domain  as  the 
UTM  domain  (though  the  true  UTM  domain  is  a  rotated 
version  of  the  PCS  coordinate  system);  we  also  call  the  two 

variables  (x,jp)  as  the  Northing  and  Easting  coordinates; 

this  terminology  is  used  to  identify  the  UTM  domain. 

The  slant-plane  where  the  image  is  formed  is  denoted 

with  (x5 ,  ys )  where  slant-plane  range, 


Slant-plane  cross-range, 

ys  =  -x  sin  0f  +y  cos  0f  (2) 

Slant-plane  ground  range, 

xg^  =  x  cos  0£+y  sin  0P  (3) 
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Target  area  mean  elevation  is  defined  by  Ztarget  and 
synthesized  flight  path  motion  angle  is  Q,  . 

We  should  point  out  that  the  2D  slant-plane  {xs,ys) 

varies  with  each  subaperture  and  should  be  denoted  with 

(  fa)  fa)  \ 

\XS  ,ys  ).  For  notational  simplicity,  we  drop  the 

superscript  l  to  identify  the  2D  slant-plane.  The  wavefront 
processor  produces  two  SAR  images  on  the  slant-plane  using 
monostatic  and  bistatic  data  for  each  Subaperture.  We  denote 

these  via  /™i0  fa, ,  y, )  and  fb(-Jt{xs,ys) .  Furthermore,  the 

coherent  signal  subspace  processing  of  these  two  imagery 
results  in  the  GSSD  that  is  the  MTI  information;  we  define 

this  by  /mtU-wJ- 

Using  the  inverse  of  the  transformations  in  equations  (1)- 
(3),  these  slant-plane  images  are  transformed  into  the  ground 
plane  of  the  PCS  coordinate  system;  we  refer  these  as 

/molds’ ys)  =>  gmonofa-fa 

/mti  CwJ  =>  s mti  h’  y ) 


III.  Residual  Clutter  Signature  Suppression 


IV.  Moving  Target  signal  model  and  image 

SIGNATURE  WITHIN  A  SUB  APERTURE 


Our  next  task  is  to  examine  the  issues  that  are  associated 
with  tracking  a  detected  moving  target  and  estimating  its 
motion  parameters.  To  do  so,  we  first  develop  a  signal  model 
for  the  signature  of  a  moving  target  and  its  image  signature 
within  a  subaperture.  In  our  model,  we  assume  the  velocity 
of  a  moving  target  is  constant  within  a  single  subaperture. 
Consider  a  moving  target  that  is  located  at  the  coordinates 


(  ^target  ’  ^target  ’  ^target  ) 


at  the  slow-time  T  —  0 ,  that  is,  at  the  midpoint  of  the 
subaperture.  Furthermore,  the  3D  velocity  of  this  moving 
target  is  denoted  with 

(  ^xtarget  ’  ^target  ’  ^ztarget  ) 

Thus,  the  radial  distance  of  the  target  from  the  radar  as  a 
function  of  the  slow-time  is 

\o  nl/2 


(t)  = 

Tadar- target  V  / 


(V 

fa 


A 

xra 


+  VW  T—  X  —V  T 

radar  '  ^xradar  L  target  ^xtarget 1 ) 


+ 


+  T  —  Y  —V  T 

radar  jradar  target  jtarget  / 

(zW  +vW  T-Z  -v 

radar  '  UradaA  U 


+ 


target  .ztarget 


(4) 


The  coherent  signal  subspace  processing  yields  an  MTI 

information  base  fy[j\  )  in  which  the  background 

s  s 

clutter  is  mostly  suppressed.  This  slant-plane  image, 
however,  does  contain  some  residual  clutter  signature.  Since 
the  slant-planes  vary  with  each  subaperture,  the  residual 
clutter  signatures  do  not  appear  at  the  same  coordinates  in  the 
subaperture  images  (though  they  correspond  to  the  effects  of 
stationary  targets).  Meanwhile,  when  the  subaperture  MTI 
images  are  transformed  into  the  UTM  domain,  that 

is, gu\i{x,y)  the  transformed  residual  clutter  signatures  do 
co-register  in  the  UTM  coordinates. 

This  property  can  be  exploited  to  further  reduce  the 
clutter  signature  in  the  UTM  domain  MTI  imagery.  For  this 
purpose,  we  construct  a  reference  image  via  a  subaperture  t 
domain  median  filtering  of  the  MTI  images  for  every  point  in 
the  UTM  coordinates.  We  denote  the  resultant  by 

<§ref  ix>y)  •  Tlhs  image  is  basically  a  representative  of  the 

common  residual  clutter  signature  in  all  the  MTI  imagery;  it 
contains  almost  no  moving  target  signature.  Next,  we  use 
this  reference  image  in  the  adaptive  filtering  method  of  signal 
subspace  processing  to  reduce  the  residual  clutter  signature  in 

each  subaperture  MTI  image  y) .  We  denote  the 

final  MTI  image  by  g|AR  (x>  >’)  • 


We  denote  the  angular  Doppler  frequency  of  this  target  via 
the  following:  ^Lim-get  )  • 

Using  SAR  signal  theory,  the  radar-target  distance  and  the 
angular  Doppler  frequency  are  related  via  the  following: 


sm 


m 

radar-target 


dr 


(*) 

vradar-target 


v<0  -v 

Vxradar  xtarget 

(v(0  -v 

Vyradar  _ytarget 


ztarget 


)(x,M 

){z, 


|  sr _  V  _ -y  sr  1  i 

radar  '  Kxradar  L  target  Kxtarget  "  )  ' 

_i_  T V  _ y  T )-|- 

radar  '  yyr adar  1  1  target  K>’target  L  )' 

(0  _!_„(*) 


4-  T  —7  —V  T) 

radar  zradar  target  ztarget  / 


R 


(d 

radar-target 


w 


(5) 


As  we  mentioned  earlier,  the  principal  assumption  in  our 
subaperture -based  approach  is  that  the  signature  of  a  moving 
target  appears  fairly  localized  in  subaperture  images.  This 
location  in  the  polar  coordinates  of  the  slant-plane  are  the 
radar-target  radial  distance  and  the  angular  Doppler 
frequency  at  the  midpoint  of  the  subaperture  (that  is,  the 
slow-time  T  —  0 ): 

lA-»w(oUii.„w(o)l  (6) 

Thus,  the  slant-plane  coordinates  of  the  moving  target  in  the 
subaperture  images  are 


xw 

Vtarget 


y 


M 

^target 


=  ^Ltarget(°)  COS^lr-,arget(0) 
=  ^radar-target (0)  S’n  41r-targe,(0) 


978- 1-4244-581 3-4/1 0/$26.00  ©2010  IEEE 


000454 


Using  the  mapping  from  the  slant-plane  to  the  UTM 
domain,  the  slant-plane  coordinates  of  the  moving  target  can 
be  converted  to  its  ground  UTM  (PCS)  coordinates: 

[^target  ?  J7 ^target  J  [^target  ?  J7 target  J 

The  distance  of  this  signature  from  the  location  of  the 
platform  at  the  mid-point  of  the  apertures  in  the  UTM  domain 
is: 


V.  Results 


The  overall  signal  processing  used  in  our  approach  is  shown  in 
Figures  1  and  2.  Figure  3  is  the  UTM  domain  reconstruction 
of  the  interrogated  for  Subaperture  1.  Figure  4  shows  the  MTI 
image  for  the  same  scene  using  the  coherent  GSSP. 


Figure  1 .  Wavefront  Reconstruction  and  Coherent  Clutter  Suppression  with 
Monostatic  and  Bistatic  Data  of  Along-Track  Monopulse  SAR 


Figure  2:  Slant-Plane  to  UTM  (Ground-Plane)  Transformation  of  MTI 
Image,  and  Residual  Clutter  Suppression 
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Figure  3:  UTM  Reconstruction:  FP-128;  Subaperture  1 
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Figure  4:  Coherent  MTI  via  GSSP  in  UTM  Domain: 
FP-128;  Subaperture  1 


Figure  5  is  the  reference  image  that  is  constructed  from  the 
median  filtering  of  the  UTM  domain  MTI-GSSD  imagery;  see 
the  block  diagram  in  Figure  2.  Figure  6  is  the  resultant  MTI 
image  for  Subaperture  1  using  non-coherent  LSSP  processing 
of  the  MTI-GSSD  image  of  Figure  6  and  the  reference  image 
of  Figure  5.  We  have  generated  movies  of  these  results  for  a 
set  of  the  subaperture. 
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Figure  6:  MTI  After  Residual  Clutter  Suppression 
FP-128;  Subaperture  1 


Figure  7a:  Detected  Moving  Targets  in  Reconstruction  and 
MTI.  FP-128;  Subaperture  1 


Figure  7b:  Detected  Moving  Targets  in  Reconstruction  and 
MTI.  FP-128;  Subaperture  1  (zoomed  version) 
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Figure  8a:  Radial  Locations  of  Detected  Moving  Targets 
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results  presented  here  were  generated  using  Gotcha  radar 
data.  This  research  will  continue  to  detect,  geolocate  and 
tracking  multiple  moving  targets  in  SAR  imagery  in  real¬ 
time. 
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